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' Abstract 

o : 

We discuss the way in which factorisation is partially maintained but nevertheless 
\ modified by process-dependent colour factors in hadronic single-spin asymmetries. 

\ We also examine QCD evolution of the twist-three gluonic-pole strength defining an 

effective T-odd Sivers function in the large- x limit, where evolution of the T-even 
e-j . transverse-spin DIS structure function is known to be multiplicative. 

: 1 Preamble 

> : 

q ; 1.1 Motivation 

^ ' . . . 

Single-spin asymmetries (SSA's) have long been something of an enigma in high-energy 

t— i ■ hadronic physics. Prior to the first experimental studies, hadronic SSA's were predicted 

to be very small for a variety of reasons. Experimentally, however, they turn out to be 

large (up to the order of 50% and more) in many hadronic processes. It was also long held 

that such asymmetries should eventually vanish with growing energy and/or p T . Again, 

however, the SSA's so far observed show no signs of high-energy suppression. 
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1.2 SSA Basics 

Typically, SSA's reflect spin-momenta correlations of the form s ■ (p A fc), where s is some 
particle polarisation vector, while p and k are initial/final particle/jet momenta. A simple 
example might be: p the beam direction, s the target polarisation (transverse therefore 
with respect to p) and k the final-state particle direction (necessarily then out of the p s 
plane). Polarisations involved in SSA's must usually thus be transverse (although there 
are certain special exceptions). 

It is more convenient to use an helicity basis via the transformation 

lt/i> = ^[l+>±i|->]- (1) 
A transverse-spin asymmetry then takes on the (schematic) form 

A (T I T) ~ (I I I) 2Im(+|-) f 
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The appearance of both |+) and |— ) in the numerator signals the presence of a helicity-flip 
amplitude. The precise form of the numerator implies interference between two different 
helicity amplitudes: one helicity-flip and one non-flip, with a relative phase difference (the 
imaginary phase implying naive T-odd processes) . 

Early on Kane et al. [1] realised that in the massless (or high-energy) limit and the 
Born approximation a gauge theory such as QCD cannot furnish either requirement: for a 
massless fermion, helicity is conserved and tree-diagram amplitudes are always real. This 
led to the now infamous statement [1]: "... observation of significant polarizations in the 
above reactions would contradict either QCD or its applicability.'" 

It therefore caused much surprise and interest when large asymmetries were found; 
QCD nevertheless survived! Efremov and Teryaev [2] soon discovered one way out within 
the context of perturbative QCD. Consideration of the three-parton correlators involved 
in, e.g. #2, leads to the following crucial observations: the relevant mass scale is not that 
of the current quark, but of the hadron and the pseudo-two-loop nature of the diagrams 
can generate an imaginary part in certain regions of partonic phase space [3]. 

It took some time, however, before real progress was made and the richness of the 
newly available structures was fully exploited — see [4]. Indeed, it turns out that there are 
a variety of mechanisms that can generate SSA's: 

• Transversity: this correlates hadron helicity flip to quark flip. Chirality conservation, 
however, requires another T-odd (distribution or fragmentation) function. 

• Internal quark motion: the transverse polarisation of a quark may be correlated with 
its own transverse momentum. This corresponds to the Sivers function [5] and requires 
orbital angular momentum together with soft-gluon exchange. 

• Twist-3 transverse-spin dependent three-parton correlators (c/. #2): here the pseudo 
two-loop nature provides effective spin flip (via the extra parton) and also the required 
imaginary part (via pole terms). 

The second and third mechanisms turn out to be related. 

2 Single-Spin Asymmetries 
2.1 Single-Hadron Production 

As a consequence of the multiplicity of underlying mechanisms, there are various types 
of distribution and fragmentation functions that can be active in generating SSA's (even 
competing in the same process): 

• higher-twist distribution and fragmentation functions, 

• /cT-dependent distribution and fragmentation functions, 

• interference fragmentation functions, 

• higher-spin functions, e.g. vector-meson fragmentation functions. 

Consider then hadron production with one initial-state, transversely polarised hadron: 

A\P A )+B(P B )^h(P h )+X, (3) 

where hadron A is transversely polarised while B is not. The unpolarised (or spinless) 
hadron h is produced at large transverse momentum PhT and PQCD is thus applicable. 
Typically, A and B are protons while h may be a pion or kaon etc. 
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Figure 1. Factorisation in single-hadron production with a transversely polarised hadron. 



The following SSA may then be measured: 

A h = da(S T ) -da{-S T ) 
T da(S T ) + da(-S T ) ' [ ' 

Assuming standard factorisation to hold, the differential cross-section for such a process 
may be written formally as (cf. Fig. 1) 

dcr = &<* ^ Xa "> f' b ( Xb "> ® d ^«'77' ® V h/c( z )i ( 5 ) 

abc aa'77' 

where f a (fb) is the density of parton type a (6) inside hadron A (B), p a aal is the spin 
density matrix for parton a, T>JJ c is the fragmentation matrix for parton c into the final 
hadron h and d<7 aa / 77 ' is the partonic cross-section: 

where M. a ^ is the amplitude for the hard partonic process, see Fig. 2. 

The off-diagonal elements of T>JJ c vanish for an unpolarised produced hadron; i.e., 
^V/c ^ ^rr' ■ Helicity conservation then implies a = a' , so that there can be no dependence 
on the spin of hadron A and all SSA's must vanish. To avoid such a conclusion, either 
intrinsic quark transverse motion, or higher-twist effects must be invoked. 



M a 




Figure 2. The hard partonic amplitude, a (3^5 are Dirac indices. 
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2.2 Intrinsic Transverse Motion 

Quark intrinsic transverse motion can generate SSA's in three essentially different ways 
(all necessarily T-odd effects): 

1. k T in hadron A requires f a (x a ) to be replaced by V a (x a , fc r ), which may then depend 
on the spin of A (distribution level); 

2. kt in hadron h allows T>JJ c to be non-diagonal (fragmentation level); 

3. k' T in hadron B requires f b {x b ) to be replaced by V b (x b ,k' T ) — the spin of b in the 
unpolarised B may then couple to the spin of a (distribution level). 

The three corresponding mechanisms are: 1. the Sivers effect [5]; 2. the Collins effect [6]; 
3. an effect studied by Boer [7] in Drell-Yan processes. Note that all such intrinsic-fc T , 
-k t or -k' T effects are T-odd; i.e., they require ISI or FSI. Note too that when transverse 
parton motion is included, the QCD factorisation theorem is not completely proven, but 
see [8]. 

Assuming factorisation to be valid, the cross-section is 

E h = ^2 [ dXa dxb d2kT d2fc ^ d ' ' 



abc act' PP'-yy' ' 



irz 



x V a (x a , k T ) p a a , a V b (x b , k' T ) p% (^) Vtf e (z, «r), (7) 



where again 

/da 
Vd£ 



) ^E^^^W ( 8 ) 

The Sivers effect relies on T-odd /^-dependent distribution functions and predicts an SSA 
of the form 

= \S T \ £ [ dx a dx b ^A T f a (x a , k%) f b (x b ) dHXa [^ kr) D h/C (z), (9) 

TTZ dt 



abc 



where Aq / (related to f^) is a T-odd distribution. 



2.3 Higher Twist 

Efremov and Teryaev [2] showed that in QCD non- vanishing SSA's can also be obtained by 
invoking higher twist and the so-called gluonic poles in diagrams involving qqg correlators. 
Such asymmetries were later evaluated in the context of QCD factorisation by Qiu and 
Sterman, who studied both direct-photon production [4] and hadron production [9]. This 
program has been extended by Kanazawa and Koike [10] to the chirally-odd contributions. 
The various possibilities are: 

d<T = ^ { G F( X a, Va) ® fb(x h ) ® d<7 ® D h/c ( Z ) 
abc 

+ A T / a (a; a ) <g> E b F (x b , y b ) <g> da' <g> D h/c (z) 

+ A T f a (x a ) <g> f b (x b ) ® da" <g> D® c {z)} . (10) 



4 



The first term is the chirally-even three-parton correlator pole mechanism, as proposed 
in [2] and studied in [4, 9]; the second contains transversity and is the chirally-odd con- 
tribution analysed in [10]; and the third also contains transversity but requires a twist-3 
fragmentation function Df~J . 



2.4 Phenomenology 

Anselmino et al. [11] have compared data with various models inspired by the previous 
possible (fc^-dependent) mechanisms and find good descriptions although they were not 
able to differentiate between contributions. The calculations by Qiu and Sterman [4] 
(based on three-parton correlators) also compare well but are rather complex. However, 
the twist-3 correlators (as in g 2 ) obey constraining relations with fey-dependent densities 
and also exhibit a novel factorisation property, to which we now turn. 



2.5 Pole Factorisation 



Efremov and Teryaev [2] noticed that the twist-3 diagrams involving three-parton corre- 
lators can supply the necessary imaginary part via a pole term; spin-flip is implicit (and 

due to the gluon). The standard propagator prescription (— •— in 
Fig. 3, with momentum fe), 
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k 2 ± is 



leads to an imaginary contribution for k 2 — > 0. A gluon with 
momentum x g p inserted into an (initial or final) external line p' 
sets k = p' — x g p and thus as x g — > we have that k 2 — > 0. 
The gluon vertex may then be factored out together with the 
quark propagator and pole, see Fig. 4. Such factorisation can be 
performed systematically for all poles (gluon and fermion): on 
all external legs with all insertions [12]. The structures are still 
complex: for a given correlator there are many insertions, leading to different signs and 
momentum dependence. 

The colour structures of the various diagrams (with the different types of soft inser- 
tions) are also different (we shall examine this question shortly). In all cases (exam- 
ined) it turns out that just one diagram dominates in the large-iV c limit, see Fig. 5. All 
other insertions into external (on-shell) legs are relatively suppressed by l/N c 2 . This 



Figure 3. Example of 
a propagator pole in a 
three-parton diagram. 
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Figure 4. An example of pole factorisation: p is the incoming proton momentum, p' the outgoing 
hadron and £ is the gluon polarisation vector (lying in the transverse plane). 
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has been examined in detail by Ramilli (Insubria U. Masters the- 
sis [13]): the leading diagrams provide a good approximation. 
The analysis has yet to be repeated for all the other twist-3 con- 
tributions (e.g. also in fragmentation). 

A question immediately arises: could there be any direct rela- 
tionship between the twist-3 and /^-dependent mechanisms? It 
might be hoped that, via the equations of motion etc., unique 
predictions for single-spin azimuthal asymmetries could be ob- 
tained by linking the (e.g. Sivers- or Collins-like) &>r-dependent 
mechanisms to the (Efremov-Teryaev) higher-twist three-parton 
mechanisms. An early attempt was made by Ma and Wang [14] 
for the Drell-Yan process, but the predictions were found not to 
be unique. Ji et al. [15] have since also examined the relation- 
ship between the /^-dependent and higher-twist mechanisms by matching the two in an 
intermediate kr region of common validity. 



Figure 5. Example of 
a dominant propagator 
pole diagram. 



3 More on Multiparton Correlators 
3.1 Colour Modification 

In [16] we provided an a posteriori proof of the relation between twist-3 and fcr-dependence. 
The starting point is a factorised formula for the Sivers function: 



dAa ~ J d 2 k T dx f s (x, k T ) e psPk ? Tr [ 7p H(xP, k T )] . 



:i2) 



Expanding the subprocess coefficient function H in powers of k T and keeping the first 
non- vanishing term leads to 



d 2 k T dx f s (x, k T ) k% e psPk ? Tr 



7p 



O H(xP, k 7 
~dkf~ 



(13) 



By exploiting various identities and the fact that there are other momenta involved, this 
can be rearranged into the following form: 



dAa ~ M dx (x) e asPn Tr 



dH(xP, k T ) 

dig, 



J k T =0 



where 



2M2' 



(14) 



(15) 



The final expression coincides with the master formula of Koike and Tanaka [17] for twist-3 
gluonic poles in high-p T processes. The Sivers function can thus be identified with the 
gluonic-pole strength T(x, x) multiplied by a process-dependent colour factor. 
The sign of the Sivers function depends on which of ISI or FSI is relevant: 



/i 1) (x) = X)C i ^ ? T(x,x) 

i 



(16) 



6 



where Cj is a relative colour factor defined with respect to an Abelian subprocess. Emis- 
sion of an extra hard gluon is needed to generate high p T and, according to the process 
under consideration, FSI may occur before or after this emission, leading to different 
colour factors. In this sense, factorisation is broken in SIDIS, albeit in a simple and 
accountable manner. Figure 6 depicts the application of this relation to high-p T SIDIS. 




Figure 6. Twist-3 SIDIS tt production via quark and gluon fragmentation. 



3.2 Asymptotic Behaviour 

The relation between gluonic poles (e.g. the Sivers function, and T-even transverse-spin 
effects, e.g. g 2 [18-22]) still remains unclear. Although there are model-based estimates 
and approximate sum rules, the compatibility of general twist-3 evolution with dedicated 
studies of gluonic-pole evolution ([23, 24] and at NLO [25]) is still unproven. 

In the large- a; limit the evolution equations for g 2 diagonalise in the double-moment 
arguments [26]. For the Sivers function and gluonic poles, this is the important kinematical 
region [4]. The gluonic-pole strength T(x), corresponds to a specific matrix element [4]. 
It is also the residue of a general qqg vector correlator by(xi,x 2 ) [27]: 

b v (xi,x 2 ) = h regular part, (17) 

x 1 - x 2 

which is defined as 

b V (xi,X 2 ) = y giAiCxi-xa)-^^ cA-Pi™^, s |^( ) ^ ^(Ax) V(A 2 )|pi, a> - (18) 

There is though one other correlator, projected onto an axial Dirac matrix: 

M*i,s 2 ) = ^J ^l^M^)^ ^ s \mil b s-D^^X^s) . (19) 

This last is required to complete the description of transverse-spin effects, in both SSA's 
and g 2 . The two correlators have opposite symmetry properties for x\ <-> x 2 : 

b A (xi,x 2 ) = b A {x 2 ,x 1 ), b v (xi,x 2 ) = -b v (x 2 ,x 1 ), (20) 

determined by T invariance. In both DIS and SSA's just one combination appears [19]: 

b-(xi,x 2 ) = b A (x 2 ,xi) - b v (x 1 ,x 2 ). (21) 
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The QCD evolution equations [20-22] are best expressed in terms of another quantity, 
which is determined by matrix elements of the gluon field strength: 



Y(x 1 ,x 2 ) = (xi - x 2 )b-(xi,x 2 ). (22) 
It should be safe to assume that b-(xi,x 2 ) has no double pole and thus 

T{x) = Y{x,x). (23) 

Evolution is easier studied in Mellin-moment form, implying double moments for 

Y(x,y): 

Y mn = J dxdyx m y n Y(x,y), (24) 

where the allowed regions are \x\, \y\ and \x — y\ < 1 (recall that negative values indicate 
antiquark distributions). We wish to examine the behaviour for x and y both close to 
unity and therefore close to each other. The gluonic pole thus provides the dominant 
contribution: 

\xmY(x,y)=Tm*) + 0(x-y). (25) 

X,J/-t-l 

In this approximation (now large m, n) the leading-order evolution equations simplify: 

ds V 2 / 

where the evolution variable is s = (3q 1 lnlnQ 2 . In terms of T{x) this is 

f W ^4(c F+ %)£ rfz [;:g ( j i)+ r W , ,27) 

which is similar to the unpolarised case, but differs by a colour factor (Cf + Ca/2) and a 
softening factor (1 — z)/(l — x). 

The extra piece in the colour factor (Ca/2) vis-a-vis the unpolarised case (just Cf) 
reflects the presence of a third active parton — the gluon. That is, the pole structure of 
three-parton kernels is identical, but the effective colour charge of the extra gluon is Ca/2. 
The softening factor is inessential to the asymptotic solution, it merely implies standard 
evolution for the function f(x) = (1 — x)T(x). For an initial f(x,Ql) = (1 — x) a , the 
asymptotic solution [28] is the same but modified by a — > a(s), with 

a(s) =a + 4:(c F + ^yjs. (28) 

For T(x), a shifts to a — 1; the evolution modification is identical; the spin-averaged 
asymptotic solutions are thus also valid for T(x). This large- a; limit of the evolution 
agrees, barring the colour factor itself, with studies of gluonic-pole evolution [23-25]. We 
should note here that Braun et al. [29] have found errors in [23] and [25]; our results for 
the logarithmic term are, however, unaffected. 
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4 Summary and Conclusions 



Viewing the Sivers function as an effective twist-3 gluonic-pole contribution [16], it is seen 
to be process dependent: besides a sign (ISI vs. FSI), there is a process-dependent colour 
factor. This factor is determined by the colour charge of the initial and final partons. 
It generates the sign difference between SIDIS and Drell-Yan at low p T , but in hadronic 
reactions at high p T it is more complicated. Such a picture is complementary to the 
matching in the region of common validity. The matching between various p T regions 
now takes the form of a p T -dependent colour factor. It also lends some justification to 
the possibility of global Sivers function fits [30]. 

We have shown that generic twist-3 evolution is applicable to the Sivers function. Its 
effective nature allows us to relate the evolution of T-odd (Sivers function) and T-even 
(gluonic pole) quantities. A vital ingredient here is the large- a; approximation, where 
gluonic-poles dominate and the evolution simplifies. The Sivers-function evolution is then 
multiplicative and described by a colour-factor modified twist-2 spin-averaged kernel [31]. 
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